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Abstract 


Continuity of local time for Brownian motion ranks among the most notable 
mathematical results in the theory of stochastic processes. This article addresses 
its implications from the point of view of applications. In particular an extension of 
previous results on an explicit role of continuity of (natural) local time is obtained for 
applications to recent classes of problems in physics, biology and finance involving 
discontinuities in a dispersion coefficient. The main theorem and its corollary provide 
physical principles that relate macro scale continuity of deterministic quantities to 
micro scale continuity of the (stochastic) local time. 


1 Introduction 


Quoting from the backcover of the intriguing recent monograph BarndorfF-Nielsen 


and Shiriaev (2010|: 


Random change of time is key to understanding the nature of various 
stochastic processes and gives rise to interesting mathematical results and 
insights of importance for the modelling and interpretation of empirically 
observed dynamic processes 


This point could hardly have been more aptly made with regard to the perspective of 
the present paper. 

The focus of this paper is on identifying the way in which continuity/discontinuity 
properties of certain local times of a diffusive Markov process inform interfacial dis¬ 
continuities in large scale concentration, diffusion coefficients, and transmission rates. 
For these purposes one may ignore processes with drift, and focus on discontinuities 
in diffusion rates and/or specific rate coefficients. This builds on related work of 
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the authors where the focus was on interfacial effects on other functionals such as 
occupation time and first passage times; see Ramirez et al. |2013 ) for a recent survey. 
We also restrict our study to one-dimensional processes, some extensions to higher 
dimensions can be found in 


Ramirez (2011). 


Dispersion problems in the physical sciences are often described by a second order 
linear parabolic equation in divergence form that results from a conservation law 
applied to the concentration of some substance. A particular class of interest for this 
paper is that of the one-dimensional Fokker-Plank equation (Q with discontinuous 
parameters D and rj and specified discontinuities in the concentration. More precisely, 
we consider the solution u{t,y) to the following problem for ?/ G M \ I, t ^ 0, 


du d Z ^ du 
dt dy dy 


op 

dy 


= 0, I3^u{t,x"r) = p u{t,x-), j€Z. (1) 


with a prescribed initial condition u{0,y) = UQ{y), y G M, and under the following 
assumptions. 


Assumptions 1. We consider a discrete set of “interfaces" 


I := {xj : j eZ}, Xo:=0, Xj < Xj+i, J G Z, (2) 

with no accumulation points. The functions D and rj exhibit jump discontinuities 
only at points in I: 

[up. := D{x+) - D{xj) PO, [rip. ■.= T]{x+)-pixj) PO, j e Z. (3) 

We further assume that D and p are functions of bounded variation in M with p 
being continuous, and D differentiable in {xj,Xj+i) for all j G Z. Finally, there exist 
constants 0 < k < K < oo such that 


k p D{x) P K, k p p{x) P K, X G M. 


(4) 


The constants {Pf} are strictly positive and such that 

> <oo. 

^ B ■ 
jez W 


(5) 


Equation (Q is often referred to as a continuity equation for the conserved quantity 
p{y)u{t,y)\ Fourier's flux law for heat conduction and the corresponding Fick's law 
for diffusion being among the most notable such occurrences. The right-side of the 
pde is the divergence of the diffusive flux \D^ of u, and the first bracket condition 
is continuity of flux at the interfaces. 

In Jean Perrin's historic determination of Avogadro's number Na from observa¬ 
tions of individual particle paths submerged in a homogeneous medium, the procedure 
was clear (Perrin 1909D . Einstein had provided a twofold characterization of the dif¬ 
fusion coefficient D in Fick's law: first as a function of variables at the molecular 
scale, including Na', and second, as the rate of growth of the variance of the posi¬ 
tion of particles in time. This meant that D, and therefore Na could be statistically 
estimated. 

If one regards (Q as a given physical law that embodies certain interfacial discon¬ 
tinuities at points in I, then the question we want to address is what corresponding 
features should be specified about the paths of the stochastic process? Our basic 
















goal is to show that the answer resides in suitably interpreted continuity properties of 
local time. 

In an informal conversation where we posed this question to David Aldous, his 
reply was that he wouldn't use a pde to model the phenomena in the first place! 
Of course this perspective makes the question all the more relevant to probabilistic 
modeling. The mathematical (probabilistic) tools are clearly available to do this, and 
much of the objective of this paper is to identify the most suitable way in which 
to express the stochastic model in relation to the underlying phenomena. It will be 
shown that the interfacial conditions at the pde level, can be characterized in terms 
of the continuity properties of a certain local time process of the associated diffusion. 
In this regard, the continuity of local time of standard Brownian motion will be seen 
to indirectly play a pivotal role. 

The evolution problem (Q can be viewed as the forward equation ^ = C*u for 
the operator C* : Dom(£*) —)■ L‘^{dy) given by 


C*f 


V dy V ^y) ’ 


for functions / G Dom(£*) satisfying, besides other decay conditions, that 


( 6 ) 


Af'AAP-f'ixj), jez. (7) 

Due to the presence of the coefficient r]{y), taking the adjoint in L^(dy) of £* 
does not generally yield an operator £ that generates a positive contraction semigroup 
on C'o(M). In fact, integration by parts yields. 




dx 


dx \r] 


and any g G Dom(£) will satisfy the following interfacial conditions: 


= 0 , 


Dix^ 

A 


( 4 ) = 


Ax] 

A 


(x. 


( 8 ) 


(9) 


We refer to the corresponding evolution problem ^ = Cv as the backwards equation. 

While physical laws are often formulated on the basis of conservation principles, 
not all such models are apriori conceived in conservation form. In fact some may 
be explicitly formulated as a specification of coefficients via a Kolmogorov backward 
equation with an operator of the form ( [lo) ) below, or directly in the form of a stochastic 
differential equation; for a variety of examples of both in the present context, see 


hydrolgy. 

Cantrell and Cosner 

(12004 1 

McKenzie et al. 

(20091; pvaskainen (|2004 

1; 

Ramirez i 

CN 

O 

CN 

Schultz and Crone (2005| biology and ecology. 

Nilsen and Sayit 


(20111 finance, Guo et al. (20101 astrophysics, and Matano and Palma (20081 from 
physical oceanography. 

To accommodate the broad class of such possible examples, the present paper 
follows the following general setting. Let £ be as in ([^ and define the operator A by 

dy' 


. 1 ^ Id 

Ag := -£(r/y) = - — 
g g dx 


dx 


The interfacial conditions satisfied by a function g G Dom(^) follow from 
can generally be written in the form: 


( 10 ) 


and 


[ 9 ]xj = 0 , Xjg\xj) = (1 - Xj)g'{xj ), j G Z 


( 11 ) 










































































for some values A := {Xj : j G Z} C (0,1) for which we further assume the following 
decay condition (equivalent to (|^ under Q): 

( 12 ) 

j&Z J 

In Section |^we construct a Feller process X = {X(t) : t ^ 0} with generator 
(Al, Dom(Al)). The significance of this association is the following: the fundamental 
solution p{t,x,y) to the backwards evolution problem = Av is precisely the 
transition probability density of X{t). Namely, for an initial condition ?;(0, x) = vo{x), 
the solution v{t,x) can be written as 


v{t,x)= [ p{t,x,y)vo{y)dy = E[vo{X{t))\X{0) 
Jr 


= x\ 


(13) 


It follows, in turn, that the fundamental solution to the original forward problem Q 
is given by 


, , r/(x) , , 

Q{t,x,y) = ■^^p{t,x,y). 


(14) 


This defines the interpretation we will use with respect to the physics of ([^ as a 
Fokker-Plank equation, and to relate it with the diffusion process X. 

Remark 1. One might note that this is the form of Doob's h-transform under the 
further constraint that y is harmonic with respect to A. However this latter condition 
is generally not appropriate for the physical examples of interest in the present paper. 
On the other hand, were it applicable it could provide an altervative approach to our 
problem via conditioning, e.g. see Perkowski and Rul\ (2012^ . 

In order to identify how the interfacial conditions ( [II| ) affect the sample paths of 
the process X, we will look at the behavior of the natural local time of X at points 
in I. This notion of local time was introduced in Appuhamillage et al. (2013) as the 
density of the occupation time operator with respect to Lebesgue measure. Namely 
for A G i3(]R), i^{t,x) is a previsible process, increasing with respect to t, such that 


lyl(X(s))ds= [ £^{t,x)dx. 

Ja 


(15) 


Our main result is the following: 


Theorem 1. Suppose D, y satisfy Assumptions I] Let X be the Feller process with 
infinitesimal generator {A, Dom(Al)) defined by (I^ [11]) and A satisfying (12). Then 


i^{t,xj) y{xJ)D{x,) A,- 


j G Z. 


(16) 


Xj ) y{xj ) D{xj ) 1 - Aj ’ 

with probability one, for any t such that x^) > 0. 

Corollary 1. Suppose D,y and j G Z} satisfy /Assumpt/ons Let u be the 
solution to the forward equation Q and X its associated Feller process. Then 

L^{t,xj) _ v{xp/3~ _ y{x+)u{t,x+) 

£^{t,xj) V{xj)/3^ y{xj)u{t,xj)' 

with probability one, for any t such that xil~) > 0. 


j G Z. 


(17) 



















The main principles to be taken from these mathematical results are as follows. 
First, Theorem and, in particular, its corollary demonstrate how the continuity 
properties of local time are reflected in the specifications of Q at interfaces. It is 
noteworthy that under the continuity of flux condition, the diffusion coefficient plays 
no role with regard to continuity of local time. In particular, its determination would 
continue to be by statistical considerations of (local) variances along the lines used 
by Perrin, while the jumps in the natural local time are a manifestation of other 
characteristics of the model. That is, the relative values of rj for given values of 
are reflected in the local time behaviors of sample paths at 
points. An example is furnished below in which the /3 and r] parameters 
are relative manifestations of geometries of both the medium and the dispersing 
individuals. 

Under continuity of flux, the continuity of local time is equivalent to continuity of 
the conserved quantity rju. In particular if ry = 1, or more generally is continuous, the 
continuity of u is a manifestation of the continuity of local time. These connections 
between continuities at the macro and micro scale are dependent on the continuity 
of flux in defining the physical model. 

A second principle arises for those contexts in which (10) is a prescribed back¬ 


, oi' vice-versa, 
interfacia 


ward equation with ry = /3 = 1, e.g. as in financial mathematics, BarndorfF-Nielsen 
and Shiriaev (2010). From a physical perspective the interface condition is not a 


continuity of flux condition, however continuity of local time occurs if and only if 


V - ^7 

D-+D 


j = l,2, 


1.1 Example: piecewise constant coefficients. 

The scope and interest of our result may be illustrated with a single interface ex¬ 
ample motivated by applications to ecological dispersion in heterogeneous media. 
Consider a population of erratically moving individuals occupying an infinitely long, 
two-dimensional duct as depicted in Figure]^ Let A{y) be the cross-sectional area 
of the duct, and 1/Ty(y) the biomass of any individual occupying the cross-section at 
a distance y from the interface ?/ = 0. Let c{t,{y,y)) denote the concentration of 
biomass, which is assumed continuous throughout, y denoting the transversal spa¬ 
tial variable. Then r]{y)c{t, {y,y)) is the concentration of individuals, for which we 
assume the following modification to Fick's law: the flux of individuals is given by 
—Ill 2 x 2 Vc, namely proportional to the gradient of the concentration of biomass. If 
D,A and ry are taken to be piece-wise constant: 


D{y) := 



if y > 0 
if y ^ 0, ’ 



if y > 0 
if y ^ 0 ’ 


then the concentration of biomass per unit length 



if y > 0 
if y ^ 0 ’ 
(18) 


u{t,y) 



(y,y))dy 


satisfies the following one-dimensional forward problem 


du 1 d f du\ 
dt 7] dy\^ 5y / 


dy 


= 0 , 


J 0 


1 

A^ 


u(t, 0 ) = 


(19) 


^n(T0+), (20) 










which is of the form Q with = IjA^. The corresponding backwards operator A 
is given by ( [lo) ) and ( [II| ) with Aq = i X denote the diffusion process 

generated by A. 


Given an initial biomass distribution uq, by virtue of ( [I^ , the solution to (20) 
can be written in terms of the transition probabilities of X as: 


u{t,y) = 


viy) 


uoix)r]{x)p{t,x,y)dx. 


( 21 ) 


Namely, p{t,x,y) is the density function of the location X{t) of individuals that 
started at x, and the paths of X may be regarded as a model for the random movement 
of the individuals. 

As in all cases of an operator A of the form (10) with piece-wise constant co¬ 


efficients at a single interface, the associated process X may be expressed explicitly 


by re-scaling an appropriate skew Brownian motion (see Appuhamillage et al. ( 2013|| 
2011); Ramirez et al. (2013); Ouknine| ( 1991D ). Specifically, 




where is skew Brownian motion with transmission probability 

A^^yr]+D~^ 


a{D,A,rj) = 


A+ y^rj+D+ -I- A- 


( 22 ) 


(23) 


The effect of the interface on the sample paths of the process X can be readily 
observed in the stochastic differential equation solved by X\ 


X{t) = 



D{X{r)) 

ri{X{r)) 


dB{r) -I- 


A+D+ - A-D- 
' 2A+D+ 


Lf(f,x), 


(24) 


where B denotes standard Brownian motion, and {L^{t,x) ; f ^ 0} is the right 


semimartingale local time at x of the unknown process X. See Revuz and Yor (1999) 
for details. For the current example, the right continuous version of natural local time 
is given by f^{t,x) := x), and for every t > 0, is discontinuous 

at X = 0 with probability one, with a jump characterized by (16). 


The nominal effect that a single interface exerts on the particle paths can now be 
elucidated by looking into further properties of skew Brownian motion, and natural 
local time. On one hand, ( [14] ) and the value of a{D,X,'n)i inform us that for an 
initial condition concentrated at y = 0, tto = 5o(dy) in (|^, the individuals will 
asymmetrically distribute on either side of the interface: 


I 71 

q{t,0,y)dy = P(X(t) > 0 X(0) = 0) = -!—a{D,X,v), 

' TjX 


for all t > 0. (25) 


On the other hand, natural local time £^{t,x) can be related to the time the 
process spends in a small vicinity of x. Of particular interest is the relative times 
particles spend at either side of zero. It follows from Corollary (|^ that: 

£l(o,,)(X(s))ds 


y+A"' 

"^0 /o 11(_,,0)(^('S)) ds V~A- 


lim 


f > 0, a.s. 


(26) 
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Figure 1: Schematics of an example of a problem leading to a one-dimensional diffusion 
process with one single interface at ?/ = 0. 


1.2 Notation and outline 

The analytical treatment here revolves around functions / : M —)> M that are measur¬ 
able with respect to the Borel cr-algebra The main function space is C'fe(M), the 

space of real valued continuous bounded functions. i?(M) denotes, in turn, bounded 
measurable functions on M. For a measure ^ on (M, i3(]R)), denotes the Hilbert 

space {/ : M ^ /r( x)^(dx) < oo}. The right and left limits of a function 
at X are denoted by f{x~^) and f{x~), and their difference is the jump operator 
[f]x = Whenever defined, the derivative of / is ^ or /', while 

f±{x) = f'{A) denote its left and right derivatives at x. 

For the general theory of one-dimensional diffusion processes used here, we refer 


the reader to Revuz and Yor (1999). A diffusion process consists of the measurable 
space (fljT") and a family of probability measures {P^; : x G M}. A sample path of 
the process is X = {Ar(t) : t ^ 0} G fl. Under P^;, the paths of X “start at x", 
namely, Pa;(X(0) = x) = 1 for all x G M. Only diffusions on M, with infinite lifetime 
are considered here. Much of the analysis of such processes is undertaken in terms 
of their scale and speed measures, which we treat as follows. For a G M, the hitting 
time of a by X is = inf{t > 0 : X{t) = a} and the scale measure s(dx) is 
characterized by 


P,(W,^ < H^) = 


X\ 


s{{a,x)) 


a < X < b. 


(27) 


s{a,b) ’ 

Every scale measure in this paper will be absolutely continuous with respect to 
Lebesgue measure, s(dx) = s'(x)dx. We also define, without room for confusion, 
the scale function ^ 

s : M —)> M, s(x) := / s'{y)dy. (28) 

Jo 

The speed measure m(dx) of X is the unique Radon measure such that 


A H^)= [ 
Jm 


[s(x Ay)- g(a)] [s{b) - s{x V y)] 
s{b) — s{a) 


m{ dy), a < X <b. (29) 


The infinitesimal generator of the process can be written in terms of the speed and 
scale measures as follows 


Yl/=-^-^/, / G Dom(Yl), 


(30) 
















in the sense that, 


-^( 2 : 2 ) -= / Af{x)m{dy) 


(31) 


Moreover, 


f{Xt) - fiXo) - / AfiXis)) ds, / G Dom(^), x G 


(32) 


is a martingale. For a given operator (,A,Dom(^)) if (32) holds, we say that X 
solves the martingale problem for A. 

The rest of the paper is organized as follows. In the next section we provide a 


construction of the diffusion process X associated to the operator (10) and identify a 
stochastic differential equation for which X is the unique strong solution. In section 
1^ we define three related but different notions of local time, including the natural 
local time and characterize its spatial continuity properties. The proof of the main 
results [I] and [I] follows directly from such characterization. 


2 On the Diffusion X 


As illustrated in Ramirez et al. (2013), the derivation of an associated process to an 
evolution operator ([^ can be achieved in several ways, including the general theory 
of Dirichlet forms, or by martingale methods. Here, we “read off” the speed and 


scale measures from the backward operator written in the form (30), and construct 


the appropriate process via a stochastic differential equation. A similar approach was 
carried out in the case of piecewise constant coefficients by Ramirez ( 2011D , and will 
be extended here to the present framework. 

Recall Assumptions on /, D and rj, and let (^, Dom(^)) be as in ([i§, 
Recursively define a sequence ipj,j G Z by 


ipj 


-_ - X = ^ j g z, ipQ := 1. 

Pj-i D{x-)Xj /37’^ 


(33) 


Then the generator (^, Dom(^)) given by (10) may be equivalently expressed by 

^ ^ ixj,Xj+i), j G Z, (34) 


r]{x) dx \ 2(pj dx j 

acting on functions in Cft(M) that are twice continuously differentiable within each 
(xj,Xj_|_i) and such that 


D{Xj ) D{x+ 




-fix.) = 


(pj 


-fixj), JGZ. 


(35) 


In view of (30) and (34), we take m(dx) := m'(x)dx and s(dx) = s'(x)dx 


with densities prescribed on M \ I by 


®^(®) ~ rui\ ’ xn'{x) = X G (xj,Xj+i), j G Z. 


D{x) ’ 


Pj 


(36) 


The existence of an associated diffusion process is established in the following 
theorem. 


























Theorem 2. Suppose D and rj satisfy Assumptions^ and A satisfies (12|. Let m 
and s be measures with densities given by ( [3^ . Then there exist a Feller diffusion 
X = {Pa; : X G M}) with speed and scale measures m and s, respectively and 

whose transition probability function p{t, x, dy) is the fundamental solution to the 


backwards equation ^ = Av with (^,Dom(^)) given by (34|, (35). Moreover, q 


in 


([I4]) is the fundamental solution to the forward problem ([^ 


Proof Note first that the boundedness assumptions on D and q, together with the 
fact that I has no accumulation points, make m and s Radon measures. Since m and 
s are assumed to have piecewise continuous densities, (30) takes the form Af{x) = 




m'{x+) ~Sx \1TT) ) Dom(^) being comprised of all functions / G C'fe(M) such that 
f'{x)/s'{x) is continuous on all of M and differentiable in M\/. Specializing to points 


in I, this specification is equivalent to the conditions in (35). The range Ran(^) of ^ 
is contained in B{R), and Dom(^) xRan(^) is a linear subset of C;,(M) x R(M). The 
existence of a diffusion process X that solves the martingale problem for (^, Dom(A)) 


could now be established under very general conditions (see for example Stroock and 


Varadhan 1979D . For our purposes however. Theorem below explicitly constructs a 
diffusion X with speed and scale measures given by m and s, which therefore provides 
a solution to the martingale problem for (^, Dom(A)) for any x G M (Revuz and 


Yor 1999 Theorem 3.12, p 308). Moreover, conditions (|^ make the boundaries 
±00 inaccessible for the process X, and it follows from [Mandl (1968 p. 38) that 
the the transition probabilities p{t,x, dy) of X make Ttf{x) = J^f{y)p{t,x, dy) a 
strongly continuous semigroup with the Feller property, namely Tt : C'fe(M) C'ft(M) 
for all t ^ 0. Let (^ 0 ; Doiii(^o)) bs the closure of the infinitesimal generator of 


{Tt : t ^ 0}. By Ethier and Kurtz] (2009 Theorem 4.1, p 182) X is generated by ^0 
in C'b(M) and is the unique solution to the martingale problem for A. It follows that 
^0 = in Dom(^o) = if F Dom(^) : Af G C'fc(M)} C Dom(^). Moreover, from 
standard semigroup theory, = Af for all / G Dom(^o). namely, p{t,x, dy) is 

the fundamental solution to the backwards equation ^ = Av. The result now follows 
from the uniqueness of fundamental solutions for parabolic differential equations (see 
for example Friedman 2013) □ 


We turn now to the construction of the diffusion X with speed and scale measures 
given by m and s in (36). The general procedure rests on the fact that the process 


s(X) is on natural scale and can be written as an appropriate re-scaling of a time- 
changed Brownian motion (see Mandl ( 1968|| ) for the general theory, or Ramirez 
(2011) for the case of piecewise constant coefficients). Then, we derive the stochastic 


differential equation that the process X should satisfy and verify that in fact, a strong 
solution exists. 

We first establish a useful lemma regarding the processes X and s(2L). 

Lemma 1. Let X be a diffusion process with scale and speed measures s and m 
defined through (|2^-(29), that are absolutely continuous with respect to Lebesgue 


measure, and have densities s' and m' respectively. 


1. Denote Y{t) = s{X{t)), t ^ 0 where s is the scale function defined in (28). 
Then Y is a diffusion with scale function sy(x) = x and speed measure my 
with density satisfying: 


m' {x) = s' {x)mY{s{x)) a.e. 


( 37 ) 











































2. The quadratic variation of X is given by: 


{X)t = 


I, m'{X{s))s'{X{s)) 


ds 


(38) 


Proof That sy{x) = x follow because Y is on natural scale (Revuz and Yor 1999 
p. 302). Let a <b and recall definition ([2^ of the speed measure. Denoting 


/ h ^ (x Ay - a){b-xV y) 

g{a, b, X, y) := - - -, x,y € (a, b). (39) 

b — a 


we can write 


E^{H^AH^)= / g{s{a),s{b),s{x),s{y))m'{y)dy. 


(40) 


Since X —>■ s{x) is an increasing function, Fx{X{t) G (a,b)} = Ps(a;)(y(f) G 
(s(a),s(6))) and the expected exit time in (40) can also be written as 


E, 


A 




s(b) 


rsib) 

ls(a) 

rb 


g{s{a),s{b), s{x),y)mY{y) dy 


(41) 


= / g{s{a),s{b),s{x),s{z))s'{z)mY{s{z))dz. (42) 


The uniqueness of the measure m (Revuz and Yor 1999 Theorem 3.6, p. 304) 
implies (37). 

To prove the second assertion, let B be Brownian motion. It follows from 
( Breiman] 1992[ Theorem 16.51) that a version of Y = s(X) can be written as 
a time change of B as follows: let (^(t) := ^mY{B{r)) dr and T{t) = 

then Y{t) = B{T{t)), t > 0. In particular, {Y)t = T{t). The quadratic variation of 
X = s“^(y) is therefore 


{X)t= [\{s-^nY{r))]^d{Y), 
Jo 


1 


/o WiXirW 


dr(r) 


By (37), and performing a change of variables, we can also write T as 


rnt) 

T{t) = 

Jo 


1 


m'Y{B{r)) 


Jq m'[X{r)) 


Combining equations (43) and (44) yields (38). 


(43) 


(44) 


□ 


The following is an extension of results of Ouknine (1991) for the case of piecewise 
constant coefficient and a single interface, and is an equation of the general type 
considered by Le Gall| ([l984[) and, more recently, Bass and Chen (2005). 


Theorem 3. Under Assumptions\^ the process X constructed in Theorem\^is the 
pathwise unique strong solution to 



D{X{s)) 

r]{X{s)) 


dB{s) 



D'iXjs)) 

2r,{X{s)) 


ds + Y, 


2A, - 1 

2A, 


L^(Txj) 


(45) 


where L^(t, x) is right semimartingale local time of X and the functions D, g are 
taken to be left continuous at points in I. 








































Proof. By Lemma [^the continuous martingale Y{t) = s{X{t)) has absolutely con¬ 
tinuous quadratic variation {Y)t = fg Z^(r) dr where Z(r) := follows 


from Karatzas and Shre^(1991 Theorem 3.4.2) that there exists a filtered proba¬ 


bility space with a Brownian motion B, such that Y{t) = fg Z{r)dB{r). Since D is 
assumed to be of bounded variation, the function can be written as the difference 
of two convex functions and 


, —In/// 1 N S''(S ^(x)) V—>. 

(s ) (dx) =--^- 77 — dx -L > 


1 


1 


s'(x 7 ) s'(x - ) 


(46) 


Applying the Ito-Tanaka and occupation times formulas (Revuz and Yor 1999 The¬ 
orem 1.5 and Exercise 1.23, Chapter VI) on X{t) = s~^(Y{t)) yields 


X{t) = 



s'_ {X {r))m'_ {X{r)) 


di3(r) — 


»"(X(r)) 


+ E 

jez L 


0 [s'{X{rWm'{X{r)) 
s'(x|) 


dr 


1 - 


s'(x. ) 


(47) 


L+{t,Xj) 


which coincides with (|4^. The pathwise uniqueness of strong solutions follows from 
Le Gall| ( |l984 ); Bass and Chen (20051, by noting that under the current assumptions, 
y/D/r] is a function of bounded variation, bounded away from zero, and the measure 
(s“^)"(dx) in (471 is finite with ^ | for all j G Z. 


□ 


Having obtained the diffusion X corresponding to the conservation form of (Q, 
in the next section we explore the role of continuity of flux in the structure of X and 
its local time. 


3 Various Notions of Local Time 


Local time has a striking mathematical role in the development of the modern theory 
of stochastic processes, from Brownian motion and diffusion, to continuous semi¬ 


martingale calculus, e.g., see Revuz and Yor 


Rogers and Williams (20001; Ito 


and MacKean (1974|. In the course of this development two particular variations on 


the notion of local time have occured as follows: 


Definition 1. Let X be a continuous semimartingale with quadratic variation {X)t. 
The right, left, symmetric semimartingale local time (smit) of X is a stochastic 
process, respectively denoted L^{t,x), L^{t,x), {t,x), x G M, t ^ 0, continuous 

in t and determined by either being right-continuous in x, left continuous in x, or by 
averaging L^{t, x) = (T+(i, x) -L L^{t, x))/ 2 , and such that in any case 


ip{X{s))d{X)s = I ip{x)L^^„{t,x)dx 


(48) 


almost surely for any positve Bore! measurable function p. 

Remark 2. The notation here is slightly different from that of Revuz and Yor (1999). 









































Observe that by choosing y? as an indicator function of an interval [x, x + e), e > 0, 
one has by right-continuity, for example, that 


L^(t, x) = lim - 
40 e 


(X(s))d(X),. 


(49) 


Similarly, L5 can be obtained by using the indicator on the interval (x — e, x]. 

The next definition is that of diffusion local time (dit) and requires Feller's no¬ 
tions of speed measure m(dx) and scale function s (x) (see and ([^). It is 
customary to define dIt only for diffusions on natural scale; e.g., Ito and MacKean 


(1974|; Rogers and Williams] (2000). Flowever, since any diffusion in natural scale is 
a time change of Brownian motion, it follows that its local times will therefore be 
themselves time changes of the local time of Brownian motion, and therefore always 
(spatially) continuous. On the other hand, for a general Feller diffusion X with scale 
function s, the transformation Y = s(X) produces a diffusion on natural scale. This 
transformation renders local time continuity as a generic property that does not fur¬ 
ther inform more specific structure of the diffusion X. We thus extend the definition 
to Feller diffusions with any scale. 

Definition 2. Let Y be a diffusion with speed measure my(dy). Then the diffusion 
local time (dit) ofY, denoted {y,t) is specified by 


T{y{r))<^r = J^ip{y)L^{y,t)mY{dy), 

almost surely for any positive Bore! measurable function ip. 

By Lebesgue's differentiation theorem, it follows that 

1 


(50) 




(51) 


with the corresponding formula for L^, and as the average. 

For the case of piecewise constant coefficients at a single interface, a particular 


local time continuity property at the interface was identified in Appuhamillage et al. 
(2013|. It was useful there to consider a modification of the more standard notions 


of semimartingale and diffusion local time to one referred to as natural local time. 
This was achieved there by exploiting explicit connections with skew Brownian mo¬ 
tion indicated above. An extension to piecewise continuous coefficients and multiple 
interfaces is obtained in the present note. 

The following modification of the definition of local time will be seen as useful in 
precisely calibrating jump size relative to the interface parameters. Just as in the case 
of semimartingale local time, one may consider right, left, and symmetric versions. 


Definition 3. Let X be a regular diffusion. The natural local time (nit) of X, right, 
left, and symmetric, respectively, denoted is specified by the occupation 

time formula 

pt p 


(p{X{s))ds= I ip{y)i^^,,{t,x)dx, 


(52) 


0 


for any positive Bore! measureable functions ip. The right and left versions are defined 
by the respective right-continuous, left-continuous versions, while the symmetric nit 
is defined by the arithmetic average of these two. 















Remark 3. In its simplest terms, the modification to natural local time is made 
physically natural by examination of its units, namely [^], whereas those of smit 
are those of (spatial) length [L], while dit is dimensionless, ffowever, as previously 
noted, its essential feature resides in the implications of continuity properties relative 
to conservation laws. In particular, this puts a notion of stochastic local time on 
par with fundamentally important principles of concentration flux and conservation 
of mass for pdes. 


The relationship between the three notions of local time above is summarized 
in the following proposition, the proof of which follows from a direct application of 
definitions and 

Proposition 1. Let X be a Feller diffusion with speed measure m and quadratic 
variation {X)t = /g g(X(s)) ds, then 

= - 1 ’ = m'{x)L^{t,x), a.s. (53) 

q{x) 


with right and left versions obtained by considering the right and left continuous 
versions of q and m' respectively. 


The celebrated theorem of Trotter et al. (1958) on the (joint) continuity of local 
time for Brownian motion is well-known for the depth it provides to the analysis of 
Brownian motion paths. This result is also naturally at the heart of the following 
general characterization of continuity of natural local time for regular diffusions. 


Theorem 4. Let X be a Feller diffusion on M with absolutely continuous speed 
measure m{ dx) = m'{x) dx and scale function s(x). Then the ratio £^{t, x)/m'{x) 
is continuous. Moreover, the natural local time of X is continuous at x if and only if 
m' is continuous at x. 


Proof. Let Y{t) = s{X{t)) = B{T{t)) as in the proof of lemma Then L^j t, x) = 
L^{T{t),x), t ^ 0, X G M. On the other hand, (53) together with ( [3^ and (Revuz 
and Yor[ |1999[ Excercise Vl.1.23), give 


i^{t, x) = 


m'{x)s'{x) X 


{t, x) = 


m' {x) 


L {t,s{x)) = 


m'{x) 


(Tit), six)) (54) 


and the assertion follows from continuity of and the scale function. 


□ 


It is noteworthy that, in general, the semimartingale local time of X is not made 
continuous by division by m' as in Theorem]^ while, as has been previously noted, the 
diffusion local time of the process transformed to natural scale is always continuous. 
From the point of view of applications the theorem shows that natural local time 
furnishes a microscopic probe to detect interfacial parameters q, when /3 = 1, or /3 
when q = D, respectively, in 0. through location and size of its discontinuities. 

The desired Theorem on the role of the continuity of flux for the process X 
defined by (34) now follows as a corollary. 


Proof of Theorem\^ . Let Xj G I. By Theorem |4| and the definitions of 
give 


m 


ifj 


in 


^^(f,x+) m'(x+) qix'^)ipj-i q{x^)D{x-)Xj 


e^it,x. ) m'(x . ) q{x- )(pj q{x- )D{x^){l - Xj) 


(55) 


□ 
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